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NOMENCLATURE 
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= accelerometer gain, Eq. (5) 

= accelerometer output,' Eq. (5) 

= moments of inertia about principal axes 
- = least common multiple of 

= torque, Eqs. (1) , (3) 

- torque, Eqs. (2), (3) 

- spacecraft spin period 

= spacecraft nutation period 

= -1/t + iQ, complex energy ■ dissipation , . 

Eq. (14) 

. = angular position of control thruster with 
respect to spacecraft axes. 

A thrust pulse duration . expressed in 
radians or degrees 

= thrust pulse duration, in seconds 

= accelerometer angular position, Eq. (5) 

= nutation phase angle at initiation of 
thrust pulse 


a 

T ■ 

03 (t) 

0 ) . 
s , 


= measure of moment of inertia ratio 
between transverse and spin axes 

= energy dissipation time constant 

‘=niufafi'oh^ ahgl?,““Eq^. “ ("4)“, “{'6)““ “ — - 

= transverse angular velocity 
= initial condition on W in Eqs. (9), (ll) /and (l4) 
= spacecraft spin rate 

= nutation frequency ‘ 
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. ABSTRACT 

This paper analytically derives the relations which determine 
the nutation effects induced in a spinning spacecraft by periodic 
precession thrust pulses. By utilizing the idea that nutation 
need only be observed just before each precession thrust pulse, 
a difficult continuous-time derivation is replaced by a simple 
discrete-time derivation using z-transforms . The analytic results 
obtained are used to develop two types of modulated precession 
control laws which use the precession maneuver to concurrently 
control nutation. Results are illustrated by digital simulation 
of an actual spacecraft configuration. 
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I. INTRODUCTION 

, Spin . stabilization is used to inertially fix a spacecraft's 
thrust vector during a transfer ellipse; however, spin about a,n 
axis other than that of maximum inertia results in an unstable 
equilibrium condition. Energy dissipation resulting from fluid 
motion in heat pipes and fuel tanks causes the nutation angle to 
increase. Furthermore, during the transfer ellipse, if it is 
necessary to process the. spacecraft's spin axis to reorient the 
thrust vector, this maneuver can also cause the nutation to increase. 

This paper analytically derives the relations which determine 
the nutation effects induced in a spinning spacecraft by periodic 
precession thrust pulses. By utilizing the idea that nutation need 
only be observed just before each precession thrust pulse,’ a 
tedious, unmanageable continuous-time derivation reduces to a simple 
discrete-time derivation which can be easily solved through the use 
of z -trans forms . . , 

The analytic results obtained are used to develop two types of 
modulated precession control laws' which use the precession maneuver 

r 

to concurrently control nutation. 

Section II briefly reviews nutation' in spinning spacecraft and 
active nutation control. Section III develops the relations which 
determine the nutation effects induced by periodic precession thrust 
pulses. Section IV develops the relations required to modulate the 
precession thrust pulse sequence so that nutation is reduced. Sec- 
tion V describes two types of modulated precession nutation control 
laws which are demonstrated in the digital simulation of a Synchronous 
Meteorological Satellite (SMS) configuration. Finally, Section VI 
presents conclusions. 


V 





II. NUTATION AND ACTIVE NUTATION CONTROL 

In 1968 Grasshoff [1] published the original -work on a 
control law which can be used to sense and remove nutation auto- 
matically. Grasshoff 's basic differential equation for the trans 
verse angular velocity of the spinning spacecraft is 

w (t) + ifio)(t) = N (1) 


and the solution for a control thrust duration of (t -t .) seconds 

is i 

is given by 

iX 






[1-e ] 


( 2 ) 


where N iri (1) is defined by Grasshoff as 


N = Nq e^^ (3) 

and X is the angular position of the thruster in the spacecraft's 
coordinate system. Thus, if a thrust pulse is initiated at the 
proper time and for the proper duration, w(t ) is less than u)(t ) 
and the transverse angular velocity is reduced. 

Grasshoff proceeds to define the nutation angle as 


0 = w(t)/(o)^-fi), (4) 

“"Nutation -is sensed by" an““accelerorrieter-w±th 'sensitive’' axis'" • 

parallel to the spacecraft's spin axis and' is removed by' properly 
timing the firing of a small thruster whose thrust axis is also 
parallel to the spin axis . . 

The sensed acceleration is given by 





(5) 


a = sin (fit + p) 

in which 0 is the nutation angle, u is the phase angle that 
defines the position of the accelerometer v/ith respect to the 
position of the thruster, and a^ is a gain that is a function 
of the accelerometer's radial position and the spacecraft's spin 
rate . 

Energy dissipation due to liquid movement in fuel tanks and 
heat pipes of the spinning spacecraft causes nutation to contin- 
ually increase. For small pertubations , the increase in nutation 
angle, 0, caused by energy dissipation is generally represented 
[3] as I 

0(t) = eQe^/"", (6) 

in which x is the time constant of energy dissipation. 

Due to the relationship in (4) between 0 and w(t), energy 
dissipation can be included in (1) to give 

u(t) + (-^ + ifi)co(t) = N (7) 

Figure 1 illustrates a ^simple, nonoptimal , nutation control 
law investigated by Taylor [2] for the International Ultraviolet 
Explorer (lEU) spacecraft. A threshold is established which rep- 
resents the amount of nutation that can be removed by a properly 
phased thrust pulse of duration equal to one spin period, i.e., 
2ir/Wg. Nutation is sensed by the accelerometer. When the next 
zero-crossing of the acceleration sinusoid occurs, the thruster 
fires and reduces the nutation to near zero. Illustrations la and 


lb show the nutation control for a small initial nutation, and 
Ic shows the series of repeated thrust pulses that are required 
for a large initial nutation. 

The nutation control law presently planned for the lUE 
operates in the same manner except that the thrust pulse is of 
a duration equal to one-half the nutation period, i.e., ir/fj. 

III. NUTATION DURING UNMODULATED PRECESSION 

Precession of the spin axis of a spinning spacecraft can be 
effected by using the same thruster that is required for nutation 
control; however, in this case, precession in a uniform direction 
is the result of a train of thrust pulses which have a frequency 
equal to the spin frequency of the spacecraft. This section shows 
how precession affects the nutation of. a spacecraft with and with- 
out energy dissipation. 

Precession Without Energy Dissipation 

Equation (1) describes the change in transverse angular 

velocity, and by virtue of (6), the change in nutation, angle .without 

energy dissipation, as a result of the forcing function N. 

-^Cons ider- -a precession maneuver consrsting^df a *^thru"stihg“ sequence 

of equally spaced thrust pulses that are At wide and have a 

period T = 2Tr/w . Thus, 
s s 

N, n*T^<t<n*Tg + At, 

w(t) + ifio)(t) = ^ n=0,l,2... 

O, n»T +At<t<(n+1)T . 

I s - s 


( 8 ) 


The continuous solution for Eq. (8) is complicated'by the 
essentially infinite pulse train. For example, the Laplace 
transform of Eq. (8) is given by 


w(s) = 


0) . 

O 


N(l-e-“=), 


S + if2 , . . /n •'Te; S. 

s (s+iJ2) (1-e s ) 


(9) 


Where Wq is the initial condition on (8) . 

(See reference [4] for a discussion of the Laplace transforms of 

pulse trains.) The transient portion of Eq. (9) has poles at 
s = 0 , s = -ifi , and s = +i2n7r/T^ , n = l, 2, 3.... 

The inverse Laplace transform of Eq. (9) is extremely complex 
and tedious because the pole at s = -ifl lies on the imaginary axis 
and cannot be separated by a contour from the poles at +i2mr/T . 
Thus, a continuous time-domain solution is formidable and not worth 
the effort because a bound on Eq. (8) will suffice. 

. The solution of Eq. (8) is straightforward if it is sufficient 
to know only the result at discrete times, for example, at the dis- 
crete point, nT , just prior to the next thrust pulse. Thus, Eq. 
(8) can be changed to the first-order, difference equation 


w (' 


o)[ (n-l)T^] 


-i^T ,, 
e s + 


N 

ifi 


,, ifiAt. -i^T 
(1-e ) e s . , 


( 10 ) 


(See reference [5] for a discussion of difference equations and 
z-transform theory.) 

^y^ taking the z-t:^r^ pj _ ( lOJ _an^^_then^^ta]y,ng__the^ 
z-transform, the discrete time solution is 


,,r ml ^-iP(nT ), N ,, iPAt, -i^T , 
wC n'^]=w^e ® In 


, -iPT 
1-e s 


-i!J(nT ) 
1-e s 



•4> 


The dynamics of (11) are more apparent If , constant parameters are combined 


to give 


le.te ' f [ li. e 




0)[nT^] = 


( 12 ) 


in which )a, is a constant given by 


, N Sin (l/2QAt) 

1 ~Q. Sin (1/2J2T ) ' 

s 


Equation (12) is periodic 


with period, I, where Z is the least common multiple of ' 


and is the nutation period 2tt/Q. In Eq. (13) , is ] 


bounded 


except for the singularity that occurs when T^=tt/Q, which is the 


case when the spin period equals one -half, the nutation period. 


Thus, nutation resulting from precession, in the absence of energy 


dissipation, is bounded except at one singularity. 


Precession With Energy Dissipation 


Equation (7) describes the change in transverse angular 


velocity, and consequently nutation, with energy dissipation. 


As in the case without energy dissipation, the use of the z- 


transform to provide a discrete-time solution offers the most 


'straightforward approach. Thus, the discrete- time solution of 


Eq. (7) is 


* ' “Ct (nT ) 

a)[nTg]=u) e~“^^'^s) + (N/a) (l-e~“^%"“'^s [ i_ + g— — (14) 

1-e s i.-e s 


where a = -(1 /t) + iQ 





in which (jj[nT^] is unbounded unless oi exactly equals k_, which ' 
s . o i . 

is only theoretically possible. ■ 

The conclusions are that with no energy dissipation and with 
a judicious choice of inertia' ratios, which determine the rela- 
tionship between T and T„, nutation resulting from unmodulated 
precession is bounded and periodic and can be kept at tolerable 
levels. With energy dissipation, nutation resulting from unmodu- 
lated precession is unbounded. 


IV. NUTATION DURING MODULATED PRECESSION 

I ■ . ' ' . ' ‘ ■ 

Because precession and nutation control can use the same 
thruster, precession thrust pulses can be selected that will con- 
currently reduce nutation. 

The effect of a single precession pulse on nutation can be 
demonstrated by a graphic representation of Eq. (1) shown in 
Figure 2a. In this figure, is the . duration in radios of the 
thrust pulse during each nutation cycle. Thus, for an initial u (t ) 

j 

at an angle , under the^inf luence._,of =.a^thr.ust^ of^'Yi^rt^ j’' “ ' 

radians in duration, to(t ) is the resultant of the fixed thrust 

vector plus m(t.) rotated through the angle A$. For the conditions 

■1 . . ' 

shown in Figure 2a, 0 ){t ,) is greater than u(t ) , i.e., nutation has 

.2 I 

increased. 

Figure 2b shows the locus of o)(t .) for a thrust pulse of given 

• 2- - 
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duration as a function of w(t ) and the initial angle $ . Note 

o 

1 

that nutation is reduced whenever w(t ) occurs in the hatched 


area . 




Figure 2c is essentially the same as Figure 2b except that 
circle A and circle B are scaled to show the initial angle 
The displacement of the scales on circle A and circle B results 
from the value of A$. Thus, for a thrust pulse of duration A$=90°, 
starting when the tip of the oj(t ) vector is at 4>^=0° on circle 


B, the result will be an a)(t, ) vector with the tip at 0° on circle 

it . 

A. Figure 2c also shows the range of for which the thrust pulse 
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physically, m is the ratio of • the radius of the thrust pulse 
lobe and the initial nutation. . - 

For a given precession maneuver, all variables in (l8). 

and (19) are fixed except for io(t . ) which is linearly related to the nutation' 

' ; ^ 

: angle : 0, at the beginning of a precession thrust pulse. 

To summarize, Eq . (18) shows that a relation exists which 

, can be used to modulate , or gate, the precession pulse train so 
that only precession pulses are allowed which concurrently reduce 
nutation. In fact, the nutation sensing accelerometer discussed 
in Section II and a threshold derived from Eq . (18) can be used 

to provide precession pulse gating. 

Figure 3a shows the use of the accelerometer output for time- 
optimal nutation control. The duration of the nutation removal 
thrust pulse is which is one-half the nutation period. Fur- 

. . thermore, any thrust, which has a pulse . duration of less than t\/Q 

but occurs during the negative half-cycle of the accelerometer 
output, will also reduce nutation, but suboptimally . 

Figure 3b is a block diagram of a modulated, precession, 
nutation" control . The end result is that only precession pulses 
=“ ^ that'^ concurrently ^reduce nutation “are “passed to “the ' thrus ter 

To substantiate that a nutation control law of the type 
shown in Figure 3 will control nutation over an extended period 
of time, it is necessary to determine the average and maximum 
time between modulated precession pulses. Because precession 
pulses are used for nutation control, the nutation removed by a 
precession pulse must be greater than nutation buildup between 
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precession pulses resulting from energy dissipation. As nutation 

increases, the gate width in Figure 3b increases to 180°. Thus, 

on the average 50 percent of the precession pulses will be passed 

and 50 percent will be blocked. Nutation can be controlled if the 

nutation increase in time 2T , where T is the spin period, is 

s s 

less than the nutation removed by one precession thrust pulse. 

The maximum time between gated precession thrust pulses is 
difficult to show analytically and depends on the ratio of T^ to 
T which may be irrational. Computer solution for the spacecraft 

O 

•r 

described in Table 1 shows that for a gate width of 50 percent, 
maximum precession pulse separation is 3T^. Furthermore, for the 
spacecraft of Table 1,. modulated precession pulses can maintain 
the nutation angle at less than 6 milliradians for precession pulse 
separation up to 50 to 60T . 

«. S 

V. MODULATED PRECESSION NUTATION CONTROL 

This section describes two control laws which use modulated 
precession pulses to control nutation. The corresponding figures 
show digital computer simulation results of these control laws 
applied to ..the Synchronous Meterological-- Satellite* (SMS) -cdnMgura- - - — 
tion specified in Table 1. 

Of the two precession modulation nutation control methods 
described belovr?, the first provides continuous precession modula- 
tion, and the second allows unmodulated precession until nutation 
increases to a threshold value at which time the precession is 
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modulated. Figure 4 shows the behavior for the continuous pre- 

/ 

cession modulation nutation control. Because the gate width , as 
given in (18), varies directly.' with nutation, the precession 
increase and nutation decrease will be fairly rapid initially. As 
the nutation reduces, fewer precession pulses are passed, and the 
increase in precession with respect to time reduces to the constant 
rate shown. . ■ • , 

Figure 5 shows the behavior for the nutation control law which 
allows unmodulated precession until nutation increases to a thres- 
hold value at -which time the precession is modulated. Thus, the 
precession and nutation increase until the nutation threshold is 
reached. At this time, precession modulation occurs. The slope 
of the precession increase diminishes as the result of • precession 
modulation, and the nutation, is held near the threshold level. 

VI. CONCLUSION 

This paper analytically derives the relations vjhich determine 
the nutation effects induced in a spinning spacecraft by periodic 
precession thrust pulses. By utilizing the idea that nutation 
need only. be-_-observed- just.«before-.-each precession^ihrust -pulse.,. - 
a tedious, unmanageable continuous time derivation reduces to a 
simple discrete-time derivation which can be easily solved through 
the use of z-transforros . 

The analytic results obtained’ are used to develop two types 
of modulated precession control laws which use the precession 
' maneuver to concurrently control nutation. In the first case the 


control law passes only those precession pulses that concurrently 
reduce nutation. Consequently, as shown in Figure 4 , 2000 seconds ■ 
are required to process 30° , but after the first ten seconds of 
the precession maneuver, nutation remains below 2 milliradians . 

In the second case the control law passes all precession 
pulses until nutation builds up to a threshold of 5.48 milli- 
radians. For nutation above the .threshold only nutation reducing 

precession pulses are passed. Precession is faster and 2100 seconds are 
required to process the entire l80°- 

By way of contrast, a simulation of a conventional control 

law which allows nutation control or precession, but not con- 
currently, for the SMS configuration described in' Table 1 • requires 
1990 seconds, to process 180° . The fuel' budget for the modulated 
precession control law shown in Figure 5 is 228.9 seconds of 
thruster fuel.' For the conventional type of precession and nuta- 
tion control the fuel budget is 307.9 seconds of thruster fuel. 

In conclusion,, for spacecraft such as those described by 
Grasshoff [1], and Taylor [2], for the SMS , and others emphasizing 
active on-board nutation control with an accelerometer as the 
nutation sensing element, the mechanization of the control laws 
described herein is straightforward and requires little additional 
electronic circuitry. 
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TABLE I 


SMS SPACECRAFT CONFIGURATION 

Mass Properties 

Mass =42.8 Slugs 

I = 244.4 Slug-Feet^ 

= 246.9 Slug-Feet^ ' 

I =97.3 Slug-Feet^ 

z 

a = 0.6039 

Angular Velocity 

w = 90 RPM (spin axis angular velocity) 

z 

= 5.6916 rad/sec (nutation frequency) 

Thruster Properties •• 

Thrust =5 lbs. 

Moment Arm = 3 Feet 

Duty Cycle = 1/12 seconds (45° of spin period) for Precession 
= 0.552 seconds (1/2 nutation period) for nutation 

Entergy Dissipation Tifne~Cbhs‘tarib" °" - - - - - 

T = 180 seconds. 
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FIGURE CAPTIONS 

Figure 1. (a) Nutation Control for Small Angles (Phase Plane). 

. (b) Nutation Control for Small Angles. 

(c) Nutation Control for Large Angles. 

Figure 2. (a) Graphical Interpretation of Equation (1). 

(b) Nutation Reduction Area 

(c) Nutation Reduction Area Scaled to Initial Angle ^ 

o 

Figure 3. (a) - Nutation Sensing Accelerometer Output 

(b) Modulated Precession-Nutation Control 

Figure 4 . SMS Simulation Results for Full Modulated Precession 
Nutation Control 


Figure 5. SMS Simulation Results for Full Precession/Threshold 
' Modulated Precession' Nutation Control 
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ABSTRACT 

> 

This NASA Grant Final Report is divided into two (2) 
parts. The first part summarizes the main effort of the grant — 
to investigate and develop precession maneuver control laws 
for single-spin spacecraft so that nutation is concurrently con- 
trolled. Analysis has led to the development of two (2) types 
of control laws employing precession modulation for concurrent, 
nutation control. Results have been verified through digital 
simulation of a Synchronous Meteorological Satellite' (SMS) con- 
figuration. ■ 

In the second part, an addition research effort, not 
originally considered, was undertaken to investigate the 
cause and elimination of nutation anomalies in dual-spin 
spacecraft. A literature search has been conducted and a 
dual-spin configuration has been simulated to verify that 
nutational anomalies are not predicted by the existing non- 
linear model. At the termination of the grant, the dual- 
spin research is still preliminary and no conclusions can be 
drawn as to the cause of the observed nutational anomalies in 
dual-spin spacecraft. 



Ill 


FOREWARD 

The research covered by this report is supported by NASA 
Grant NCR 26-003-069, "Nutation Control During Precession of 
a Spin-Stabilized Spacecraft". The work was performed by Dr. 

o ... 

Javin M. Taylor, Principal Investigator, Richard J. Donner , 
Master of Science Candidate in Electrical Engineering, and 
Vehbi Tasar, Master of Science Candidate in. Computer Science , 
at the University of Missouri - Rolla under the supervision 
of Dr. Thomas. W. Flatley, Earth Observations Systems & Systems 
Engineering Division, Space Applications and Technology 
Directorate, Goddard Space Flight Center. 

This report represents the Final Report for this project 
as required by NASA Provisions fpr Research Grants. 
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PART I. NUTATION CONTROL IN SINGLE-SPIN SPACECRAFT 

The purpose of Grant NGR 26-003-069 was to investigate 
variations in the basic precession maneuver control laws 
for single-spin spacecraft so that nutation is concurrently 
controlled. Nutation -induced by precession was analyzed in 
detail. Subsequently, two control laws employing precession 
modulation for concurrent nutation control were developed 
and verified through digital simulation of a Synchronous 
Meteorological Satellite (SMS) configuration. Results were 
reported in the Semi-Annual Status Report "Nutation Control 
During Precession of a Spin-Stabilized Spacecraft" in December 
1973. 

Manuscripts describing the research performed in this part 
of the grant effort have been submitted to the American in- 
stitute of Aeronautics and Astronautics Journal of Spacecraft 
and Rockets in their Synoptic/backup paper format. Appendix 
-A contains the Synoptic manuscript and Appendix B contains 
the backup paper manuscript, both entitled, "The Use of Pre- 
cession Modulation for Nutation Control in Spin-Stabilized 
Spacecraft". 

The grant to the University of Missouri - Rolla, through 
direct financial aid or research involvement, supported the 
obtainment of the Master of Science in Computer Science degree 
by Mr. Vehbi Tasar and the Master of Science in Electrical 
Engineering degree by Mr. Richard J. Donner. 
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PART II. NUTATION ANOMALIES IN DUAL-SPIN SPACECRAFT 

•As an additional research effort, not originally con- 
sidered in Grant NGR 26-003-069, a preliminary study was 
undertaken to investigate the cause and elimination of 
nutation anomalies in dual-spin spacecraft. The intent was 
to look at anomalous nutation behavior in dual-spin space- 
craft, in general, and, in particular, with regard to three 
NASA spacecraft, ITOS-D, ITOS-F, and AE-C. This behavior * 
is characterized by either the decay or growth of induced 
nutation to some non-zero terminal value, 

A literature search was undertaken and a partial bib- 
liography is included in Appendix C. The published litera- 
ture on dual-spin dynamics generally fails to provide in- 
sight into mechanisms which might be responsible for this 
anomalous nutation behavior. However, [16] describes an 
effort to determine the cause of apparently similar behavior 
for a particular dual-spin spacecraft. Mechanical cross- 
coupling due to dissipation in the bearing assembly connecting 
' . ' '' 

the two spinning bodies was identified in this "case as the 
probable cause. 

Yet another study [20] pin pointed the cross-coupling 
mechanism as electronic- in nature. Here nutation influenced 
the position sensor and produced errors in the despin control 
torque . 

Communications with NASA regarding ITOS-D, ITOS-F, and 
AE-C have also suggested the interface between the spinning 




bodies as the area wherein anomalous nutation behavior is 
generated. Prior to investigating the interface model, the 
dual-spin configuration described by Phillips [20] was 
simulated to verify that nutational anomalies were not pre- 
dicted by the existing nonlinear model. No nutational 
anomalies were observed in the simulation of the dual-spin 
model of [20] . 

At the termination of Grant NGR 26-003-069, the dual- 
spin research was still in a preliminary stage and no conclu- 
sions could be drawn as to the cause of the observed nutational 
anomalies in dual-spin spacecraft. 
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NOMENCLATURE 


I ,I 
X y z 

N 




a 

X 




= . accelerometer gain 
= accelerometer output 

= moments of inertia about principal axes 
= torque 

spacecraft spin period 

= spacecraft nutation period 

= -1/t +. in, complex energy dissipation 

= angular position of control thruster 
with respect to spacecraft axes. 

= thrust pulse duration, in angle of 
nutation cycle 
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= thrust pulse duration, in seconds 

= accelerometer angular position 

= nutation phase angle at initiation of 
thrust pulse 

= measure of moment of inertia ratio 
between transverse and spin axes 

= energy dissipation time constant 

= nutation angle- 

= transverse angular velocity 

.= spacecraft spin rate 

= nutation frequency 

THEME 


Spin stabilization is used to inertially fix a spacecraft's 
thrust vector during a transfer ellipse;, however, spin about an 
axis other than that of maximum inertia results in an unstable 
equilibrium condition. Energy dissipation resulting from fluid 
mot-ion in heat pipes and fuel tanks causes the nutation angle to 
increase. Furthermore, during the transfer ellipse, if it is 
necessary to process the spacecraft's spin axis to reorient the 
thrust vector, this maneuver can also cause nutation to increase. 

This paper analytically derives relations, which determine the 
nutation effects induced in a spinning spacecraft by periodic- 
precession thrust pulses. By utilizing the idea that nutation need 
only be observed just before each precession thrust pulse, a tedious 
continuous time, derivation reduces to a simple discrete-time derivation 
which can be easily solved through the use of z-transforms . 

The analytic results obtained are used to develop two types of 
modulated precession control laws which use the precession maneuver 



to concurrently control nutation. 


CONTENTS 


In 1968 Grasshoff [1] published the original work on a control 
law which can be used to sense and remove nutation automatically. 
Grasshoff 's basic differential equation for the transverse angular 
velocity of the spinning spacecraft is 

o)(t) + iflto(t) = (1) 

and the solution for a control thrust duration of (t,-t ) seconds 

, 1 o 

is given by 

„ i X . 

X „-ifi(t, -t„) . *^0® rn -ifi(t,-t„), 

w(tj^) = co^e 10 + — £^ --- ■ [1-e 1 0 ] (2) 

Grasshoff defines the nutation angle as the peak value of 

6 = w(t)/(w^-fi) . (3) 

•S' 

Nutation can! be sensed by an accelerometer with sensitive axis 
parallel to the spacecraft's spin axis and can be removed by 
properly timing the firing of a small thruster whose, thrust axis 
is also parallel to the spin axis. The sensed acceleration is 


given by 


a = aQ0 sin (f2t + y) . 


Energy dissipation causes nutation to continually increase. 
For. small pertubations , the increase in nutation angle, 9, is 
generally represented [3] as 

0(t) = (5) 


in which t is the time constant of energy dissipation. 
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Due to the relationship in (3) between 0 and u)(t), energy 
dissipation can be included in (1) to give 

i X - 

( 6 ) 


w(t) + (-— + ifi)u)(t) = NqG 


Precession of the spin axis of a spinning spacecraft can be 
effected by using the same thruster that is required for nutation 
control. For a precession maneuver consi’sting of a thrusting sequence 

t 

of equally spaced thrust pulses that are At wide and have a period 

T = 2u/ui , (1) becomes .. . ' 

s . s 


(jo(t) + iQoj(t) = > 


N, n*T <t<n*T + At, 
' s- s 


n=0 ,1,2... 


(7) 


0, h*T +At<t< (n+l)T. 

s - s 


The continuous solution for (7) is complicated by the 
essentially infinite pulse train. The inverse -Laplace transform 
of the Laplace transform of (7) is difficult because the pole at 
s = “in cannot be separated by a contour from the poles at ±i2mr/T^, . 
n=l,2,3,.... (See reference [4] for a discussion of Laplace trans- 
forms of pulse trains . ) 

The solution of (7) is straightforward if it is sufficient 
to know only the result at discrete times, for example, at -the dis- 
crete point, k T^ , just prior to the next thrust pulse. Thus, (7) 
can be changed to a first-order difference equation. (See reference 
[5] for a discussion of difference equations and z-transform theory.) 
By taking the z-transform and then taking the inverse z-transform. 


the discrete time solution for (1) is 

N e^^ 

. -ifi(nT 0 , , ifiAt) -ifiT . 

w[n T]=u)Qe s +^-^(l-e e s[ 


, -ifiT 
1-e s 


^-ifi(nT^) 

+ - — ^ ^ ] 

, ifiT \ 
1-e s 


( 8 ) 


The discrete time solution for (6) can be similarly obtained as 



(9) 


N G ^('1 ( n HP ^ 

1-e s 1-e s 


where a =_-(l/x) + in. 

The conclusions are that with no energy dissipation and with 
a judicious choice of inertia ratios, which determine the relation- 


ship between and nutation resulting from unmodulated' pre- 

cession is bounded and periodic and can be kept at tolerable levels. 
With energy dissipation, nutation resulting from unmodulated pre- 
cession is unbounded. Thus, some sort of nutation control is required 
during a precession maneuver. . 

Because precession and nutation control can use the same 
thruster, precession thrust pulses can be selected that will con- 
currently reduce nutation. The effect of a single precession pulse 
on nutation can be demonstrated by a graphic representation of (2) 
shown in Figure la. 

The range of the initial position <i>^-can be analytically derived 
from (2). After trigonometric manipulation, replacement of by • 

(jjo and (t^-t^) by At, the expression for the range of for 

which u) (t^) <u) (t^) is given by 


arccos(-m sin (At/2) ] (At/2 ) <tQ< 

arccos(-m sin (At/2) ] ^^^-X- (At/2) , (10) 


in which arccos [*]jj arccos refer to angles in the II 

and III quadrants, respectively, and m is given by 


m = (Nya) (l/o)(t^) ) . (11) 

To summarize, (10) shows that a relation exists which can be 
used to modulate, or gate, the precession pulse train so that only 
precession pulses are allowed which concurrently reduce nutation. 




In fact, the nutation sensing accelerometer discussed previously 
and a threshold derived from (10) can be used to provide precession 
pulse gating. 

Figure lb shows the use of the accelerometer output for time- 
optimal nutation control. Any thrust which has a pulse duration of 
less than t\/Q but occurs during'' the negative half-cycle of the 
accelerometer output, will also reduce nutation, but suboptimally . 

Figure 'Ic is a block diagram of a modulated precession, nuta- 
tion control. The end result is that only precession pulses that 
concurrently reduce nutation are passed to the thruster. 

Two precession modulation nutation control methods are now 
presented. The first provides continuous precession modulation, 
and the second allows unmodulated precession until nutation increases 
to a threshold value at which time the precession is modulated. The 
corresponding figures show digital computer simulation results of 
these control laws applied to the Synchronous Meterological Satellite 
(SMS) configuration specified as follows : mass, properties : mass = 

42.8 slugs, = 244.4, = 246.9, = 97.3, all in slug-feet^, 

a = 0.6039; angular velocity: = 90 rpm, = 5.69 rad/sec;' 

thruster properties: thrust = 5 lbs., moment arm = 3 feet, duty 

cycle for precession = 1/12 second (45° of spin period) , duty cycle 
for nutation = 0.552 seconds (1/2 nutation period); energy dissipation 
time constant: t = 180 seconds. 

Figure 2a shows the behavior for the continuous precession modu- 
lation nutation control. Because the gate width, as given in (10), 
varies inversely with nutation, the precession increase and nutation 
decrease will be fairly rapid initially. As the nutation reduces. 
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fewer precession pulses are passed, and the .increase in precession 
with respect to time reduces to the constant rate shown. Consequently, 
as shown in Figure 2a, 2000 seconds are required to process 30°, but 
after the first ten seconds of the precession maneuver nutation remains 
below 2 milliradians . 

Figure 2b shows the behavior for the nutation control law which 
allows unmodulated precession until nutation increases to a thres- 
hold value of 5 .48 'milliradians at which time the precession is modu- 
lated. Thus, the precession and nutation increase until the nutation 
threshold is reached. At this time, precession modulation occurs. 

The slope of the, precession increase diminishes as the result of 
precession modulation, and the nutation is held near the threshold 
level . 

In the second case the control law passes all precession pulses 
until nutation builds up to a threshold of 5.48 milliradians. Conse- 
quently, precession is more rapid. As figure 2b shows, 2100 seconds 
are required to process the entire 180°. 

By way of contrast, a simulation of a conventional control law 
which allows nutation control or precession, but not concurrently, 
for the SMS configuration described requires 1990 seconds to process 
180° . The fuel budget for the modulated .precession control law shown 
in Figure 2b is 228.9 seconds of thruster fuel.* For the conventional 
type of precession and nutation control the fuel budget is 307.9 seconds 
of thruster fuel. 

In conclusion, for spacecraft such as those described by Grass- 
hoff [1] , and Taylor [2], for the SMS , and others emphasizing active 
on-board nutation control with an accelerometer as the nutation sensing 
element, the mechanization of the control laws which use the precession 


maneuver to concurrently control nutation is straightforward and 
requires little additional electronic circuitry. 
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Figure 1. 


Figure 2 . 


FIGURE CAPTIONS 

(a) Graphic Representation of Equation (2) . 

(b) Nutation Sensing Accelerometer Output. 

(c) Modulated Precession Nutation, Control. 

(a) SMS Simulation For Full. Modulated Pre- 
cession Nutation Control. 

(b) SMS Simulation. For Full-Precession-to- 
Threshold Modulated Precession Nutation 
Control . 
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